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LMXB AND IMXB EVOLUTION: I. THE BINARY RADIO PULSAR PSR J1614-2230 
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ABSTRACT 

We have computed an extensive grid of binary evolution tracks to represent low- and intermediate 
mass X-ray binaries (LMXBs and IMXBs). The grid includes 42,000 models which covers 60 initial 
donor masses over the range of 1 — 4 Mq and, for each of these, 700 initial orbital periods over the 
range of 10 — 250 hours. These results can be applied to understanding LMXBs and IMXBs: those 
that evolve analogously to CVs; that form ultracompact binaries with Porb in the range of 6 — 50 
minutes; and that lead to wide orbits with giant donors. We also investigate the relic binary recycled 
radio pulsars into which these systems evolve. To evolve the donor stars in this study, we utilized a 
newly developed stellar evolution code called "MESA" that was designed, among other things, to be 
able to handle very low-mass and degenerate donors. This first application of the results is aimed at 
an understanding of the newly discovered pulsar PSR J1614-2230 which has a 1.97 Mq neutron star, 
-Porb = 8.7 days, and a companion star of 0.5 M©. We show that (i) this system is a cousin to the 
LMXB Cyg X-2; (ii) for neutron stars of canonical birth mass 1.4 Mq, the initial donor stars which 
produce the closest relatives to PSR J1614-2230 have a mass between 3.4 — 3.8 Mq; (iii) neutron stars 
as massive as 1.97 Mq are not easy to produce in spite of the initially high mass of the donor star, 
unless they were already born as relatively massive neutron stars; (iv) to successfully produce a system 
like PSR J1614-2230 requires a minimum initial neutron star mass of at least 1.6 ± 0.1 Mq, as well as 
initial donor masses and Porb of ~4.25 ± 0.10 Mq and ~49 ± 2 hrs, respectively; and (v) the current 
companion star is largely composed of CO, but should have a surface H abundance of '^10 — 15%. 
Subject headings: stars: binaries: general — stars: evolution — stars: pulsars: individual (PSR 
J1614-2230) — accretion, accretion disks — X-rays: binaries 



1. INTRODUCTION 

There has been an-ongoing effort to understand the 
evolution of low- mass X-ray binaries (LMXBs) since their 
basic nature was under stood back in the late 1960's and 
early 1970's (see, e.g. , Faulkner! 119711: iRappaport et al.l 
1983t iWebbink et all Il983l: iJoss fc RappaportI |1984|: 
Nelson et al.' '1985"; Pvlvser & Savoniic 1988, 1989'; 
Bhattacharva & van den He uycl 19 91; Iben ct ah 1995;; 
Podsiadlowski et al. 2002; P fahl et al. 2003; N elson" et al.l 
2004i : iBelczvnski et al.i.2008f ). There are two rather dis- 
tinct parts of the evolution to consider: (i) the formation 
of a neutron star in orbit with a low-mass donor star, and 
(ii) the subsequent portion of the evolution when mass 
is transferred from the companion donor star to the neu- 
tron star. 

The difficulty with the first part of the evolution in- 
volves the conceptual problem of keeping the binary sys- 
tem bound while the massive progenitor of the neu- 
tron star (NS) explodes in a supernova event. A 
part of this problem was addressed by invoking a com- 
mon envelope phase during which the lower-mass sec- 
onda ry ejects the envelope of the massive NS progen- 
itor (jPaczynskil 119761: iMever fc Mever-Hofmeisteii 119791: 
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iWebbinkI 119791: iBhattacharva fc van den Heuvell I199H) . 

However, even this solution is not so straightforward 
in that there may be insufficient gravitational energy 
release between the inspiraling low-mass star and the 
core of the NS progenitor to su ccess fully eject t he en- 
velope (see, e.g.. lDewi fc Taurisl[2000HPfahl et al.l l20031. 
This difficulty could be overcome by invoking secon- 
daries that are of intermediate mass, e.g., 2 — 4Mq 
to incr ease the gravitational energy release (see, e.g., 
iPfahl e t al. 2003). In turn, there has been a persis- 
tent conceptual misunderstanding that the transfer of 
mass from a star of 2 — 4Mq onto a 1.4 Mq NS was 
dynamically unstable. That was shown to be untrue 
fPvlvser fc Savonij ^IgSS", Igsgl; 'Tauris fc SaTO^Il99l 
Podsiadlowski fc Rappaport 2000; Ki ng et al.l 12001 : 
PRP02; see also IDavies fc Hansenlll998t (King fc Ritted 
IT99I. 

Starting in the 1970's and continuing until the present, 
there have been numerous evolution studies of a lim- 
ited number of LMXBs and IMXBs exploring the var- 
ious paths which lead to very different intermediate and 
end-stage pr oducts. One of the more syst ematic was 
the st udy bv [Podsiadlowski et al.l ()2002t ) and IPfahl erall 
(|2003f l which involved about 150 LMXB and IMXB sys- 
tems. Given the two-dimensional parameter space of ini- 
tial Porb,i and M2,i (where Porb.i and M2,i are the initial 
orbital period and mass of the donor star), the types of 
evolutionary paths for LMXBs or IMXBs is really quite 
large and varied. Possibilities include evolution to a min- 
imum orbital period of ~70 min with very low-mass H- 
rich donor stars; evolution to an ultracompact state with 
Porb in the range of 6 — 50 minutes and He-rich donors; 
and evolution to wide orbits of days-to-months with low- 
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mass giant donor stars. 

In order to better understand the many possible evo- 
lution paths of LMXBs and IMXBs in a more system- 
atic way, we have calculated an extensive grid of binary 
models encompassing 42,000 initial combinations of Porb.i 
and Af2,i- This is two orders of magnitude larger than 
the study we conducted in 2002. We took advantage 
of a newly developed stellar evolution code whose equa- 
tions of state allow for the evolution of very low mass 
stars with cold dense interiors, and a cluster of comput- 
ers which speeds up the overall calculation beyond what 
was readily available a decade ago. 

With the recent discovery of PSR J1614-2230, with 
the most massive neutron star known (1.97 ± 0.04 Mq), 
a relatively close 8.7-day orbit, and a fairly massive white 
dwarf companion (0.5 Mq), we decided to first apply our 
evolution calculations to understanding the origins of this 
system. The specific goals are to understand whether a 
natal NS with canonical mass of ~ 1.4 Mq can accrete suf- 
ficient material to grow to nearly 2 Mq , to see whether 
the donor star mass is consistent with the observed Porb, 
and to investigate how this system is related to possible 
progenitors in the guise of the LMXB - Cyg X-2. In 
all, we found ^500 of our evolution tracks which pro- 
duce systems that are at least generically related to PSR 
J1614-2230. 

In this paper we introduce our set of binary evolution 
calculations which cover an extensive grid of initial or- 
bital periods and donor masses (|J2])- In 5j3]we show how 
our binary evolution calculations apply directly to PSR 
J1614-2230. In g] we discuss a number of the general 
lessons we have gleaned from this study. 

2. OVERVIEW OF LMXB AND IMXB BINARY EVOLUTION 

2.1. Binary Evolution Calculations 

In this work we start with binary systems with an 
unevolved companion star in a circular orbit with 
an already formed neutron star of mass 1.4 Mq. 
The prior evolution, leading to the NS and low- or 
intermediate-mass donor star does not concern us in this 
study. However, this prior phase of evolution involves 
a much more massive progenitor of the NS (perhaps 
8 — ISMq), as well as a common envelope phase that 
unveils the He/CO core of the progenitor which, in 
turn, evolves to collapse and the formation of a neu- 
(see, 
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The subsequent evolution of the incipient LMXB or 
IMXB is computed with a combination of a newly de- 
veloped Henyey code called "ME SA" (Modules for E xper- 
iments in Stellar Astrophysics; IPaxton et al.l [20Tol) and 
a binary driver code. The advantages of MESA for our 
calculations are that (i) it is specifically designed to in- 
clude equations of state that are able to handle the low- 
mass and hi gh degree of degene racy that the donor stars 
achieve f see IPaxton et al1l2010l for details), and (ii) it is 
robust when used in a binary evolution code that allows 
for 'hands-off' evolution through all the various stages of 
the donor star and phases of mass transfer. We found 
that none of the evolution models failed to run to their 
expected completion. 



The binary driver code evolves the binary orbit, includ- 
ing the effects of mass transfer and angular momentum 
losses due to gravitational radiation, magnetic braking, 
and possible mass ejection from the system. It keeps 
track of where the Roche lobe is in respect to the atmo- 
sphere of the donor star being evolved by MESA, and de- 
cides how much mass to remove from the donor star dur- 
ing each evolution time step. It also manages the various 
disparate timescales involved in the donor star as well as 
in the evolution of the binary orbit (details of the driver 
code will be presented in a future work, bu t its operation 
is ver y similar to that described in Ma dhusudhan et ahl 
l2008f) . For the specific models presented here, we as- 
sumed that for sub-Eddington mass transfer rates, 90% 
of the mass transferred is retained by the neutron star, 
but that for higher transfer rates the accretion onto the 
NS was limited to that set by the Eddington limit. Mat- 
ter ejected from the system was assumed to carry away 
the specific angular momentum of the neutron star. 

No tidal evolution (between the donor star's rotation 
and the orbital angular momentum) was incorporated in 
the code. We regard this as typically a < 10% effect 
on the evolution, both before mass transfer commences 
and afterward. Magnetic braking is in cluded (accord- 
ing to eq. [36] in iRappaport et al.lll983D for all stars ex- 
cept those that are completely convective or have a mass 
> 1.4 M0. As an approximation, we also assume that 
magnetic braking operates even in those cases where the 
donor star underfills its Roche lobe and continuous syn- 
chronization between the donor star and the orbit is not 
guaranteed. We have verified by numerical experiment 
that this approximation makes little difference in the fi- 
nal distribution of evolution tracks shown in Fig.[TJ In 
this work we do not consider either thermal-ionization 
disk instabilities or the X-ray irradiation of the donor star 
during the course of the evolution calculations. Their 
effects can be approximated via the application of after- 
the-fact algorithms - which we discuss in a future work. 
However, for the purposes of the calculations emphasized 
in this paper, neither effect is very important. 

We computed binary evolution tracks for 42, 000 
LMXBs and IMXBs over a grid of 60 initial donor masses 
uniformly distributed over the range of 1 — 4 Mqj and 
700 initial orbital periods distributed over the range of 
10 — 250 hours in equal logarithmic steps. With our dense 
coverage of the initial binary parameter space, most of 
the basic types of LMXB and IMXB evolutionary tracks 
are explored. 

The "initial orbital period" in these runs is defined as 
Porb following the birth of the neutron star. The orbits 
are assumed to be circular, and the donor stars start 
on the ZAMS in our evolution calculations. The binary 
evolution is considered 'complete' when either (i) 10 Gyr 
have elapsed, or (ii) the mass transfer becomes dynami- 
cally unstable (this later condition occurs for only ~2.4% 
of our systems, typically the ones with the largest initial 
donor masses and orbital periods) . After the envelope of 
the donor star has been transferred to the neutron star 
(or ejected from the system), MESA continues to evolve 
the relic He or He/CO core of the donor until a total 
elapsed time of 10 Gyr has passed. 



2.2. Results 



LMXB and IMXB Evolution 
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Fig. 1. — Superposed set of 38,000 binary evolution tracks in the P^ih ~ ^2 plane. Each time a track crosses one of the 1300 X 1200 grid 
points in the image, the evolution time step is stored in that pixel. The color scaling is related to the logarithm of the total system time 
spent in each pixel. The accumulated evolution times (per pixel) range from 10^ yr for purple to 10^" yr for red. The initial grid of models 
is visible in the upper right portion of the plot, outlined by a red box. 




Fig. 2. — Evolution tracks of all 515 systems that terminated their evolution with Porb = 9 ± 1 day, and M2 = 0.5 ± O.IMq. The colors 
span cumulative evolution times (per pixel) of between 10^ yr (light blue) to 10^ yr (green). The small yellow region represents H-burning 
for up to 3 X 10® yr (per pixel). The location of PSR J1614-2230 is marked with a red square; that of Cyg X-2 with a red diamond. The 
four white tracks represent other typical evolutionary paths for L/IMXBs (see text for details,) 
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A graphic presentation of our 38,000 binary evolu- 
tion tracks, representing LMXBs and IMXBs is shown 
in Fig. [1] in the Porb — M2 plane. We divide this plane up 
into 1300 X 1200 discrete pixels. As the evolution tracks 
pass through this plane, we record, in a cumulative fash- 
ion, the evolution time that is spent in crossing a given 
pixel. When all the tracks have been co-added in this 
fashion, we display the logarithm of the total accumu- 
lated time in a pixel with color shading. The cumulative 
dwell times per pixel range from 10^*^ yr for the red re- 
gions to 10^ yr for the purple. The red box in the upper 
right of the diagram outlines the starting grid of initial 
models. 

To better understand what the different types of 
LMXB and IMXB evolutions contribute to diagram, we 
show in Fig. [2] a set of illustrative individual tracks, la- 
beled "CV", "UC", "Gl", and "G2". The CV-like tracks 
are analogous to the evolution of cataclysmic variables, 
where the donor star is not very evolved at the time 
when mass transfer commences, remains H-rich, and gets 
eaten away by mass transfer until the orbital period 
reaches a minimum (at ^70 min.), after which the orbit 
slowly expands (see, e . g..lPaczvns ki fc Sienkiewicz 1981i; 
iRappaport et al.lll983l : lHowell et al.ll2001L and references 
therein) . These systems commence mass transfer for Porb 
shorter than the so-called 'bifurcation period' of between 
~20 — 30 hrs (depending on the initial donor mass). 
The ultracompact systems ("UC") start mass transfer 
in somewhat wider orbits, very near the bifurcation pe- 
riod, with the donor stars correspondingly more evolved, 
typically with H just having been depleted at the stellar 



center (see. e.g.,lNel son et al."1986': ' Tutukov et al 
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1989; 



Podsiadlowski et ai.l 120021: iNelson fc RappaportI l2003f ) 



Because of the high He content in the core of the donor 
star, these systems are able to evolve to much shorter 
orbital periods (some as short as 6 min.) Finally, those 
systems that commence mass transfer in even wider or- 
bits (e.g., "Gl" and "G2"), lead to the formation of a 
well-defined He core and to an increase in orbital period. 
For those donor stars which start with A/2 < 2.2^/0, 
their degenerate cores prescribe a mass-radius relation, 
which in turn dictates an orbital period-core mass re- 
lation, Porb(-^c), which is clearly delineated by the lo- 
cus of green terminus p oints in Fig.[T] which fol l ows the 
PnrhiMr) relation (e.g., Raupaport et al. 1995; ' Ergmal 
[TOOa iTauris fc Savo niic 1999;' Nelson et al. 2004). 

For donor stars with initial masses > 2.2 Mq, how- 
ever, this Porb(-^c) relation does not have to be fol- 
lowed, and their evolution typically terminates consid- 
erably before this Porb(-Wc) boundary is reached. This 
is compounded by the fact that, due to the high mass 
ratios involved, the initial phases of the mass transfer, 
while dynamically stable, are governe d by the thermal 
timescale of the donor star (see, e.g., IKing et al.ll200ll: 
iPodsiadlowski fc RappapOTtI I2OQOI; PRP02; for further 
references see i)3.2l) . These rates can therefore become 
greater than 10~^Mo yr~\ which is far in excess of the 
Eddington limit for a N^. Therefore, the mass trans- 



® As high as these rates are, they typically remain below the 
rates required for the on set of 'hypercritical accretion' (see, e.g., 
IHouck fc ChevalierlUQQlT ) which might have allowed for the growth 
of higher mass neutron stars. 



fer is highly non-conservative, and most of the mass is 
ejected from the system. The latter class of systems, we 
believe, is responsible for the evolution of both Cyg X-2 
and PSR 1614-2230. 

In a forthcoming paper we will concentrate on various 
other aspects of these evolutionary paths and describe 
in greater detail the various types of LMXBs and recy- 
cled pulsars that they lead to. For the present work we 
focus on the evolutionary paths leading to Cyg X-2-like 
systems and PSR J1614-2230. 

3. BINARY EVOLUTION OF THE PSR J1614-3320 SYSTEM 

3.1. Properties of PSR J 1614-3320 

PSR J1614-2230 (iDemorest et al.ll2010D is a millisec- 
ond radio pulsar in an 8.7 day orbit around a (presum- 
ably) predominantly carbon-oxygen (CO) white dwarf 
(0.5±0.006 Mq). The recent measurement of the Shapiro 
delay in PSR J1614-2230 has yielded a NS mass of 
1.97 ± 0.04 Mf;^ , making it the mos t massive pulsar 
known to date (jPemorest et al.ll2010[) . Given the con- 
stituent masses and Porb — 9 days, this recycled pul- 
sar system shows all the signs of having descended from 
an LMXB very much like Cyg X- 2 (Pp^h ^ 9.8 days 
Mns - 1.8 M0 and Mo - 0.6 Mq; iCasares et £01119981 : 
lOrosz fc KuulkersI 119991 ). The orbital eccentricity of 
~10~^ indicates that the evolution involved an extensive 
period of mass transfer via a Roche-lobe filling donor 
star. The characteristic spin-down age is ~5 Gyr. 

3.2. Evolution Tracks Leading Close to PSR J 1614-2230 

Because of the very interesting and unique system 
parameters of PSR J1614-2230, we were motivated to 
search our grid of tracks for all those systems which 
terminate their evolution with Porb = 9 ± 1 day and 
M2 = 0.5 ± 0.1 Mq. In all, there were 515 such systems; 
their evolution in the Porb— M2 plane is shown in Fig. [2] as 
light blue and green tracks. Note that they all start with 
donor masses in the range of 3.35 < M2,i < 3.75 Mq, and 
Porb.i in the range of 2 — 4 days. 

Note how the light blue tracks all decrease in Porb un- 
til the donor mass reaches ~2Mq, after which Porb in- 
creases again until mass transfer has ceased when the 
donor star no longer has an extended envelope. The ini- 
tial decrease in Porb results from the fact that mass is be- 
ing transferred from the more massive donor to the less 
massive NS. The total duration of the light blue portion 
of the tracks corresponds to a characteristic evolution 
time of only ~ a Myr. This rapid phase of mass trans- 
fer is referred to as "thermal timescale" mass transfer 
since it takes place on the thermal timescale of the ra- 
diative donor star that is more massive than the accreting 
NS (see, e.g., IPodsiadlowski fc Rappaport"2000'; PRP02; 
for other related references see below in Wi.2}i . Once 
the donor masses reach ~0.8 Mq the thermal timescale 
mass transfer is over, and the evolution slows down to 
timescales of tens of Myr - the nuclear evolution timescale 
of the donors (see the green and yellow portions of the 
tracks in Fig. [2]). 

We have selected one of these 515 tracks, as illustra- 
tive, in that its end product best matches the properties 
of PSR J1614-2230. We show details of that particular 
evolution in Fig.|3l The various panels show the evolu- 
tion of Porbi -^2, P2eff, M2, M, and information on the 
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Fig. 3. — Detailed evolution vs. time of several parameters associated with the system that best represents the evolution of PSR J1614- 
2230 for a NS with initial mass 1.4 Mq (the initial time in the plots, to, was chosen as 262 Myr to emphasize the mass-transfer phase). The 
top panels show the evolution of Porbi ^2, and T^g as functions of time. The bottom panels show the evolution of stellar mass, M, and 
central composition for the donor star as functions of time. In the mass panel (d) the blue and red curves are for the donor star and NS, 
respectively. In the composition panel blue, cyan, green, and red are the mass fractions (at the center of the donor star) of H, He, C12, 
and 016, respectively. The black curve is the surface H abundance. 

chemical composition of the donor (star 2), as functions 
of the evolution time. A reference time of 262 Myr has 
been subtracted from the evolution time to enhance the 
short-timescale features around the mass-transfer phases. 
Note that M refers to the mass loss rate from the donor 
star, 

The short-lived drop in Porb (panel a) corresponds to 
the thermal timescale mass transfer event (the large spike 
in M in panel e). In panel (f) we see that the central H 
mass fraction is only ~6% by the time mass transfer com- 
mences. The relatively long plateau in M (panel e) at 
~10~^ Mq yr~^ which lasts for ~35 Myr, corresponds to 
the phase of mass transfer driven by nuclear evolution 
during which time the remainder of the H in the core 
and much of the envelope is largely consumed. Once a 
He core has formed, there is an approximately 4-Myr in- 
terval during which H-shell burning occurs. This causes 
the radius of the donor star to expand from ^3.5 Rq to 
~6i?0, during which time the mass transfer rate ranges 
between 10~^ — 1O~^M0 yr~^. The orbital period in- 
creases from '^3 days to ^^9 days, and the donor-star 
mass drops to 0.47 M©. It is during this second phase 
of mass transfer when the accreting NS undergoes its 
largest growth in mass. After the H-shell burning phase, 
the donor star contracts, Roche lobe contact and the 
mass-transfer phase is over, and the He in the core of 
the donor burns to C and O. The companion, which 
ends up being 90% C and O (with 10% of the mass 
in a He envelope), contracts to its final degenerate ra- 
dius of 0.015 i?©. This entire evolution is sometimes re- 
ferred to as "case AB" since the mass transfer occurs 
near the end of the main sequence. (For a closely related 
evolution scenario for Cyg X-2 see King fc RitteHITggl 



iPodsiadlowski &: RappaportI 120001 .) and is not the same 
as the mass accreted by the NS. 

To put the evolution of PSR J1614-2230 into some con- 
text, we show in Fig. 2] a plot of the end points of all 
the binary evolutions we ran in the plane of the final 
neutron-star mass vs. the final donor-star mass (now a 
white dwarf) - with Porb,f restricted to > 10 hr. In all, 
there are some 14,000 system end points represented in 
this figure. Even though the points form patterns and 
lie along pseudo tracks, it is important to note that they 
are not evolution tracks, but rather end points of numer- 
ous evolutions. These patterns result from the discrete 
nature of the grid of starting values for Porb.i and 
The location of PSR J1614-2230 is marked with an or- 
ange square. 

The end point locations in this Mns,f — M2,i plane are 
color coded according the final value of the orbital period. 
The heavy green dots correspond to Porb,f in the range 
of 8.7 ± 1 day, roughly commensurate with Porb of PSR 
J1614-2230. The vertical column of heavy green dots 
originates from systems with initial donor-star masses 
< 2.2 Mq which follow the core-mass radius relation (see 
discussion below), and where the mass transfer rates are 
well below the Eddington limit. Thus, the NSs in these 
systems are able to grow substantially. For higher-mass 
donors (i.e., > 2.2 Mq) that commence mass transfer 
near the end of the main sequence (case AB), the final 
NS masses lie in the range of ~2.1 — 1.65 Mq with white 
dwarf masses inversely correlated with the NS mass, and 
lying between ^^0.26 and 0.50 A/q (the continuing line of 
heavy green dots). The line of systems with the most 
massive white dwarfs (between 0.5 and 0.65 Mq) arise 
from case B mass transfer (where the initial donor masses 
were > 2.2 Mq), but their companion NS masses are all 
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Fig. 4. — The final neutron-star mass vs. the final mass of the 
white dwarf companion for different ranges of final Porbi a-s indi- 
cated by the color scaling (in hours). The dark green points are the 
systems with a final P^-ch most appropriate to match PSR 1614- 
2230. The orange square marks the location of PSR J1614-2230, 
while the large yellow circle indicates the system plotted in detail in 
Fig-El Along the locus of dark green points, the higher the initial 
donor mass the lower the final neutron star mass. The bottom- 
most line of systems is due to case B mass transfer with initial 
donor masses > 2.2 Mq. The higher mass neutron stars are very 
likely unphysical, and are a consequence of the high mass capture 
fraction assumed for sub-Eddington accretion, and not allowing 
these stars to collapse to black holes. 

< 1.5 Mq. 

This 'competition' between the NS mass and that of 
the residual white dwarf companion is due to the fact 
that, in general, the more the donor star is evolved (hence 
larger core masses) when mass transfer commences, the 
higher is the thermal-timescale mass-transfer rate, and 
therefore the lower the mass retention fraction by the NS 
(where the Eddington limit is greatly exceeded). More 
specifically, the case B mass-transfer systems have only 
a thermal-timescale mass-transfer phase, and there is no 
sub-Eddington phase during which the NS mass can grow 
significantly. Therefore, raising an initial NS mass of 
1.4 Mq to 2Mq remains a difficulty for the model (see 
also PRP02). 

The relation between the final orbital periods and the 
white dwarf companion masses of the systems ending 
with Porb > 10 hours is shown in Fig. [5] Again, these are 
14,000 end points, and not evolution tracks. The location 
of PSR J1614-2230 is marked. The systems lying along 
the leftmost boundary of these end points (marked with a 
solid line) all represent He white dwarfs of mass between 
~0.17Mq and ~0.28Mq. These systems form from 
donor stars of initial mass < 2.2 Mq, and are well rep- 
resented by a nearly unique Pnrh(Mr) relation (see, e.g., 
Rappaport et all 119951: lErgmal fl99a iTauris fc Savoniid 
1999( ). The dashed line to the left of this set of end points 
is the expression given bv lRappaport et all ()1995[ ) which 
was already known to systematically somewhat under- 
estimate the white dwarf mass for systems with short 
orbital periods below which the relationship was not de- 




0.2 0.4 0.6 O.i 

Final White Dwarf Mass (Mq) 

Fig. 5. — The final orbital period vs. the final mass of the 
white dwarf companion star for ~14,000 systems whose final or- 
bital period lies in the range of 10 — 1000 hr, or ^1/2 to 40 
days. The color coding of the dots is as follows: cyan and blue 
are for M2,i > 2.2 Mq and cases B and AB, respectively (see 
text); the green and red represent M2,i between 1.4 — 2.2 Mq 
and 1.0 - IAMq, respectively. The location of PSR J1614-2230 
is marked with an orange square. The dashed purple line is the 
short-period, low mass end of the Rappaport et al. 1995 Porb(A/c) 
relation, while the solid purple curve is a fit to the He white 
dwarfs (at the left boundary) and also subsumes the Rappaport 
et al. P^^^{Mc) relation for white dwarf masses in the range of 
0.3 — 1.2 Mq which is off this plot. This steeper slope predicting 
higher-mass w hite dwarfs for a given Porb wa-s earlier alluded to by 
I'l'auris fc B avoniic ( 1999). 

signed to work. The solid line is given by 



orb 



(1 + 25m3-5 + 29m6)3/^ 



days 



(1) 



where nic is the white dwarf mass expressed in solar 
units. We devised this expression to fit both the results 
shown in Fig.[S] for low- mass white dwarfs and to incor- 
porate the expression of iRappaport et all ()1995f ) which 
extends all the way to white dwarfs of mass 1.2 Mq. 

The systems distinctly to the right of this Porb (Mc) 
relation originate from donor stars with Af2 ^ 2.2 Mq 
and do not follow the core mass-radius relation for gi- 
ants. Note the "gap" in these more massive white dwarfs 
(which runs approximately at ~— 45°) dividing systems 
that evolved in so-called case AB (on the left side) 
from those that formed from somewhat more evolved 
donors at the time mass transfer commences, i.e., case 
B. These white dwarfs become progressively more CO- 
rich in composition towa rd the higher masses. ( For other 
related st udies see, e.g., Ilben fc Tutu kov"198 5| PRP02; 
Tauris & Savoniid 



|20 



2000; Nels on et al] 



Ilia iHan e'raiTT2000: T auris et al.l 



200^^ 



Finally, in this regard, we note that PSR J1614-2230 
is quite representative of the systems which started with 
donors of mass ^2.2 Mq and ended up with white dwarfs 
of mass in the range 0.26 — 0.65 Mq. The mass transfer 
rate in case B systems (toward the right) is so high that 
the neutron stars do not have much chance of growing 



LMXB and IMXB Evolution 



7 



to interestingly high masses. By contrast, the systems 
to the left of the 'gap' between cases AB and B have 
intervals of near- or sub-Eddington rates which allow the 
highest mass NSs to be grown (see Fig.|4|). In general, 
the higher-mass neutron stars tend to be found with the 
lower-mass white dwarfs. 

3.3. Neutron Stars with Higher Initial Masses 

As discussed in the previous section, none of our evolu- 
tion tracks, starting with neutron stars of canonical mass 
1.4 Mq, produced the observed 1.97 M© NS with the req- 
uisite combination of orbital period and white dwarf mass 
to match the PSR J1614-2230 system. Neutron stars as 
massive as 2Mq, and higher, were produced with the 
correct orbital period, but not in combination with a 
massive enough white dwarf. Similarly, the orbital pe- 
riod and white dwarf mass combination is easy to repro- 
duce, but not with the correct Porb- We find that the 
final values of M2 and Mns are anticorrelated when Porb 
is fixed to 8.7 ± 0.5 day. The difficulty with producing 
high-mass NSs in orbit with massive white dwarf com- 
panions is that the progenitors of these white dwarfs are 
initially substantially more massive than the NS, result- 
ing in very rapid, thermal-timescale mass transfer (at 
rates as high as 10~^ Mq yr~^). Unless this is followed 
by a substantial interval of sub- or near-Eddington accre- 
tion rate, the neutron star will likely be prevented from 
accreting a significant fraction of the donor star's mass. 

One obvious solution to this problem of producing 
more massive NSs is to start with higher natal mass neu- 
tron stars. To this end, we have run a series of smaller 
subgrids of 700 additional starting models with initially 
higher-mass neutron stars. We find that the minimum re- 
quired initial mass neutron star for which we could reach 
a final NS mass of 1.97 Af© turned out to be 1.6 Mq. 
The starting donor mass for these more 'successful' mod- 
els, however, increased to ^^4.25 ± 0.10 M© (compared 
to ~3.8Mo), the initial orbital periods remained near 
^^49 ± 2 hours, and the final companion CO white dwarf 
mass was ^0.49 Mq, somewhat closer to the observed 
value. An illustrative evolution sequence that produces a 
~2 Mq NS, and otherwise closely resembles PSR J1614- 
2230, is shown in Fig.[6l 

Finally, in regard to setting a limit on the mass of 
the natal neutron star, we note an important caveat. 
Throughout our calculations we adopted a very simple 
prescription for the Eddington-limit rate of accretion 
onto the neutron star of A/Edd = 3 x 10^^ Mq yr^^, inde- 
pendent of the neutron star mass or radius, or the chemi- 
cal composition of the accreted material. For most of the 
evolutionary phases we covered in our 42,000 tracks, the 
accretion rate is either below 10~^ Mq yr~^, in which 
case the value we choose for A^Edd is unimportant, or 
it so high, during the thermal timescale mass-transfer 
phases, that the choice of M-^dd also does not affect the 
very small fraction of mass that can be retained by the 
neutron star. However, for the production of systems 
that resemble PSR J1614-2230, with its massive neutron 
star and other specific properties, the largest growth of 
the neutron star occurs during the second phase of mass 
transfer in case AB evolution where M is in the critical 
range of 10~^ — 1O~^M0 yr~^, and where the value of 

A^Edd very much affects how much the neutron star can 
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Fig. 6. — Detailed evolution vs. time of several parameters asso- 
ciated with the system that matches the parameters of PSR J1614- 
2230 for a NS with initial mass 1.6 Mq and donor star of initial 
mass of 4.25 Mq. As in Fig.|3] to is an initial time that is sub- 
tracted off so as to emphasize the mass-transfer phase. The top 
panels show the evolution of the stellar masses, and the radius of 
the donor star, while the bottom panels show the evolution of M, 
and the central composition of the donor star. The color coding 
in the mass curves and chemical composition are the same as in 
Fig.[3l The final masses of the neutron star and white dwarf in this 
sequence were 1.95 Mq and 0.49 Mq, respectively, and the final 
orbital period (not shown) was 206 hours (8.4 days). 

grow. Adjusting the value of MEdd to be "more correct" 
is complicated by the following issues: (i) the unknown 
radius of the NS (by a factor of ^50%), (ii) general rela- 
tivistic corrections (of ~20%), (iii) unknown allowed fac- 
tors by which the Eddington limit may be violated in NS 
accretion (of ~ factors of 2) , and (iv) changes in chemical 
composition of the accreted material - which is especially 
important as the material becomes He-rich later in the 
evolution. Also, the amount of mass gained by the NS 
is only some ~80% of what was transferred through the 
accretion disk due to the loss of rest mass that is ra- 
diated away. Taking all these factors into account, we 
have run a range of tests and have concluded that such 
effects and their uncertainties, lead to an uncertainty in 
the minimum mass for a neutron star to reach 1.97 Mq 
of ~±0.1 Mq. Therefore, the range of required minimum 
neutron star masses could conceivably be expanded to 

1.6 ±0.1 Mq. 

4. SUMMARY AND CONCLUSIONS 

We have used MESA to compute a large grid of 42,000 
binary evolution models for LMXBs and IMXBs. We 
showed the broad sweep of possible evolutions, from sys- 
tems which attain orbital periods as short as 6 minutes 
to those which grow to long-period binaries with giant 
donor stars. We leave a detailed discussion of the bulk of 
these results for a future paper. Here we have focused on 
what we can learn about the evolutionary paths to the 
newly discovered binary pulsar, PSR J1614-2230. We 
have selected a subset of the evolution models (515 in 
total) which lead to systems like PSR J1614-2230 and its 
evolutionary cousin, Cyg X-2, to examine in more detail. 
In particular, we show how an orbital period of 8.7-days 
can be easily understood, as can the 0.5 Mq compan- 
ion star. We show that for the proposed scenario, the 
degenerate companion star is mostly (i.e., 90%) C and 
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O (with a surrounding shell of He which comprises 10% 
of the white-dwarf mass), but there is also a thin outer 
envelope that is composed of ~15% H. 

We found, however, that starting with neutron stars of 
canonical mass 1.4 Mq, we were not able to produce the 
observed 1.97 Mq NS with the requisite combination of 
orbital period and white dwarf mass to match the PSR 
J1614-2230 system. We were able to evolve neutron stars 
as massive as 2Mq, with the correct orbital period, but 
not in combination with a massive enough white dwarf. 
Also, the orbital period and white dwarf mass combina- 
tion was easy to generate, but not with the correct Porb- 
As discussed above, the difficulty with evolving high- 
mass NSs in orbit with massive white dwarf compan- 
ions is that the progenitors of these white dwarfs are ini- 
tially substantially more massive than the NS, resulting 
in very rapid, thermal-timescale mass transfer (greatly 
in excess of the Eddington limit), and the neutron star 
is thereby prevented from accreting a significant fraction 
of the donor star's mass. 

At this time, we are therefore tentatively led to con- 
clude that the initial mass of the NS would have to have 
been higher than the canonical value of 1.4 M0 in or- 
der to declare that the binary evolution of this system is 
fully understood. By running some 700 supplementary 
models with initially higher-mass NSs, we found that to 
successfully produce a system like PSR J1614-2230 re- 
quires a minimum initial neutron star mass of at least 
1.6 ± 0.1 Mq, as well as initial donor masses and Porb of 
~4.25 ± 0.10 Mq and --49 ± 2 hrs, respectively. 

For completeness in regard to producing high neutron- 
star masses, we point out that a number of the dynami- 
cally stable case B evolution tracks that we have gener- 
ated yield mass transfer rates in excess of several times 
10~^Mq yr~^ that last for intervals of several thousand 
years. If 'hypercritical accretion', where the gravitational 
energy is carried off in neutrinos and the Eddington limit 
is thereby circum vented, is able to occur during these in- 
terva ls (see, e.g.. iHouck fc ChevalieH[T99lt iBrown et al.l 
120001 ) ■ then perhaps there is a chance for the neutron 
stars to grow during this phase of the evolution. How- 
ever, it is not clear to us that the hypercritical accre- 
tion scenario, which was developed for spherical accretion 
onto a neutron star, would be applicable to accretion via 
a disk (but, see lMoreno Mendez et "alll2008[ ). 

A more general conclusion from this L/IMXB study 
is that there is a subclass of systems which start with 
intermediate mass donor stars (of > 2.2 Mq), with Porb 
well above the 'bifurcation period' (in the range of 2 — 4 
days), which leads to systems like Cyg X-2 and PSR 



J1614-2230. Such systems terminate their mass transfer 
before the donor star develops a degenerate core, and 
they end up with Pnrh well belo w the Pnrh(M^) relation 
/See alsollben fc Tutukovl [l985l: PRP02;[HaneraD[2000; 
iTauris et al.ll2000t ). 

Another important general lesson that we can take 
from this broad look at L/IMXB evolution is that 
intermediate-mass donor stars can evolve to virtually all 
the known types of LMXB systems that exist at the cur- 
rent epoch. These include CV-like evolution paths, ul- 
tracompact X-ray binaries, and systems with giant donor 
stars. This is an important finding because it is signifi- 
cantly easier for intermediate-mass stars to successfully 
eject the envelope of their massive companion progenitors 
of the NSs, and then remain bound during the ensuing 
supernova explosion. 

Finally, we have provided a perspective on where in 
the Porb and white dwarf mass plane we can expect to 
find the He white dwarfs that follow the Porb(Afc) rela- 
tion, as well as where the systems with the most massive 
neutron stars should be found. Depending on the up- 
per mass limit to a neutron star, M^s,r:na.x, Fig.|4] shows 
that there should be a fair number of black holes with 
mass between Mns^max and ~2.8 Mq in binaries with He 
or He/CO white dwarfs that range between ~0.2 and 
0.4 Mq and with orbital periods in the range of 1 — 40 
days. One way to detect such important relics of stellar 
evolution is to search for white dwarfs with interestingly 
high orbital velocities (i.e., v ~ 300 (Porb/days) km s~^) 
and unseen companions. An important caveat to this 
is that if such low-mass black holes form in significant 
numbers in LMXBs that are still undergoing mass trans- 
fer, they would be directly detected by their accretion. 
However, there is no evidence for LMXBs with low-mass 
black holes (see, e.g., iFarr et al.ll20lol: lOzel et al.l[201Q) . 
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